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ABSTRACT 


1.  Preparing  Instltutioni  RUTGERS  -  TOE  STATE  UNIVERSITY 

College  of  Engineering 
New  Brunswick,  New  Jersey 

2.  Title  of  Report!  Evaluation  and  Development  of  a  Rational 

Theory  for  the  Design  of  Sewage  Stabilization 
Ponds  t 

3.  Principal  Investigator:  Marvin  L.  Granstrom,  Ph.D. 

4.  Number  of  pages  ^  53}  Illustrations  -  16{  and 
date  -  January  30,  1963* 

5.  Contract  Number!  DA-49-193-MD-2317 

6.  Supported  by: 

U.  S.  Army  Medical  Research  and  Development  Command 
Department  of  the  Army 
Washington  25,  D.  C, 

The  purpose  of  this  study  was  to  evaluate  a  method 
suggested  for  the  design  of  sewage  oxidation  ponds.  The 
suggested  method  appeared  to  be  not  applicable  when  tested 
by  a  large  amount  of  operational  data.  This  writer  proposes 
that  oxidation  ponds  are  usually  overdesigned  and  that  beyond 
some  detention  time,  or  below  some  loading  value  termed 
’’critical  loading  value,”  the  biological  activity  of  the  pond 
is  primarily  self-'perpetuating  and  cyclic,  and  does  not  serve 
to  reduce  BOD  or  conform  bacteria.  It  is  suggsted  that 
operational  experience  be  considered  in  light  of  this  con¬ 
cept  to  lead  to  a  uniform  and  systematic  accounting  of  the 
performance  of  oxidation  ponds.  A  brief  discussion  of  the 
problem  of  odors  due  to  anaerobiasls  is  included. 
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EVALUATION  AND  DEVELOPMENT  OF  A  RATION^L  THECXIY 

FOR  THE 

DESIGN  OF  SEWAGE  STABILIZATION  PONDS 

The  work  reported  herein  wae  aade  possible  by  a  Research 
Contract  No,  DA-49-I93-MD-23I7  between  Headquarters y  U,  S,  Aray 
Medical  Research  and  Developaent  Coaaand,  Office  of  the  Surgeon 
General  and  Rutgers,  the  State  University.  The  study  was  conducted 
by  Professors  W.  H.  Schllmaeyer,  J,  P.  Lawler,  and  M.  L,  Granstroa 
of  Rutgers  University,  The  latter  served  as  Project  Director.  The 
The  period  of  the  contract  was  1  July  through  31  August  1962, 

I.  INIRCDUCTION 

In  recent  years  there  has  been  an  expanding  Interest  In  the  use 
of  a  staple  detention  basin  as  part  of,  or  as  a  ooaplete  sewage 
treataent  facility.  Such  a  facility  has  been  teraed  oxidation  nond 
or  etabillaatlon  pond.  The  foraer  will  be  used  In  this  report.  The 
detention  periods  ooaaonly  used  vary  froa  tens  to  hundreds  of  days. 

The  designs  that  have  been  used  are  priaarlly  eaplrlcal  In 
nature  but  several  atteapts  have  been  aade  In  the  recent  past  to 
suggest  appropriate  paraaeters  for  the  actual  design  of  oxidation 
ponds.  With  respect  to  structure  aost  writers  suggest  that!  the 
pond  be  enclosed  by  a  dike,  Inflowing  surface  water  be  excluded,  the 
top  of  the  dike  be  at  least  eight  feet  wide  to  allow  aaohlne  opera¬ 
tion,  the  waterside  slope  be  3  or  4  to  one,  the  dike  should  be 
planted  with  grasses  to  provide  erosion  protection  and  to  prevent 
growth  of  deep  rooted  plants  and  weeds,  the  Interior  bank  of  the 
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dike  could  be  paved  to  prevent  erosion  and  where  seepage  Is  ex¬ 
cessive  the  bottoB  and  sides  could  be  sealed,  and  the  bottom  should 
be  essentially  level  and  cleaned  of  vegetation  prior  to  putting  the 
pond  In  operation*  There  seemed  to  be  some  variation  In  suggested 
Inlet  and  outlet  structures.  The  simplest  Inlet  structure  proposed 
and  used  Is  a  horizontal  pipe  discharging  horizontally  at  a  point 
approximately  In  the  center  or  at  least  200  feet  from  any  bank  In 
the  pond.  Other  suggestions  were  that  the  Inlet  structure  be  similar 
to  those  commonly  found  In  a  settling  tank.  The  outlet  structure 
might  be  an  overflow  weir  with,  In  cases  where  level  control  Is 
desired,  and  It  usually  Is,  some  arrangement  for  selection  of  efflu¬ 
ent  at  different  depths.  Furthermore,  the  effluent  should  come  from 
a  level  several  inches  below  the  surface  of  the  water  to  prevent  ex¬ 
cessive  algae  carryover, 

Hermann  and  Gloyna^^^  have  suggested  a  formula  for  the  compu¬ 
tation  of  volume  of  a  lagoon  for  a  given  loading  as  follows! 

V  s  5.37x10“®  nqy  1,072®®“'*’ 

V  »  acre  feet 

N  =  population  equlv, 

q  =  gpcd 

y  =  5  day-20O  BOD 

T  =  operating  temp,,  °C, 

Also,  In  another  paper^^^  they  suggest  that  BOD  loadings  in 
Ibs/acre/day  in  the  climatic  regions  similar  to  Austin,  Texas  be 
based  partially  on  the  final  disposition  of  the  effluent  as  follows! 
as  land  Irrigation  water,  200;  Into  a  diluting  stream,  150;  and  Into 
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an  Intermittent  stream,  50.  They  further  suggested  thati  the  depth 
vary  from  two  to  three-and-one-half  feet,  Intra-pond  recirculation 
was  not  warranted,  ponds  should  operate  in  series  or  in  parallel, 
the  Influent  be  evenly  distributed  across  the  end  of  the  pond  and 
several  outlets  be  provided. 

Oswald,  Gotaas.  ,  Golueke  ,  and  Kellen^^^  developed  design 
equations  to  allow  determination  of  volumes  and  depths  necessary  to 
maintain  aerobic  conditions,  at  least  at  the  upper  depths  of  the 
pond,  with  appropriate  consideration  of  the  available  light  inten- 
tenslty  and  strength  of  the  waste.  They  considered  an  overdesigned 
pond  to  be  possibly  inefficient  because  an  excess  of  oxygen  pro¬ 
duction  by  algae  would  raise  the  pH  too  high  for  good  biological 
activity . 

(4) 

Neel  et.  al  have  suggested  for  ice-free  locations  that 
allowable  BOD  loadings  in  Ibs/acre/day  could  be  computed  by  dividing 
the  lowest  monthly  average  of  langleys  by  two.  A  langley  is  defined 
as  a  gram  calorle/sq.  cm.  of  incident  radiation.  They  listed  the 
low-monthly  average  of  langleys  at  a  number  of  cities  in  this  hemi¬ 
sphere  . 

A  most  significant  paper  by  TJarals  and  Shaw  entitled  "  A 
Rational  Theory  for  the  Design  of  Sewage  Stabilisation  Ponds  in 
Central  and  South  Africa”  has  been  developed  in  the  recent  past^^^. 
Their  paper  is  divided  into  five  sections.  In  section  I  the  funda¬ 
mental  differential  equation  governing  the  concentration  of  BOD  and 
faecal  bacteria  in  a  pond  is  derived.  Various  solutions  of  equations 
are  given  and  the  relationship  between  the  kinetic  activity  of  a 
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river  and  a  aerlea  of  ponda  is  established.  In  section  II  experi¬ 
mental  evidence  is  presented  to  verify  and  to  indicate  the  limita¬ 
tions  of  the  theory.  In  section  III  criterion  within  the  framework 
of  the  theory  is  developed  to  determine  the  maximum  loading  on  an 
aerobic  pond  before  anaerobic  conditions  develop.  In  section  IV 
the  theory  is  applied  for  the  development  of  a  design  procedure 
for  a  series  of  ponds.  In  section  V  brief  consideration  Is  given 
to  the  kinetics  of  recirculation  in  ponds.  It  is  believed  by  this 
writer  that  the  paper  Is  of  Importance,  The  hypotheses  and  design 
criteria  established  are  the  subjects  of  this  present  report. 

This  present  report  Is  divided  Into  four  sections.  Section  I 
is  the  introduction.  Section  II  is  a  presentation  and  discussion 
of  the  data  In  relationship  to  the  derived  equations.  The  data 
used  in  this  report  came  from  a  variety  of  sources  in  the  United 
States,  The  data  considered  In  this  section  are  the  Biochemical 
Oxygen  Demand  <B,0,D,),  9  day-20°C,t  and  conform  counts  of  sewage. 
Section  III  Includes  a  discussion  of  the  discrepancies  between  the 
hypothesized  mathematical  model  and  the  actual  data,  A  suggestion 
is  made  for  an  alternative  mathematical  model.  Section  IV  consists 
of  some  suggestions  on  areas  for  future  studies, 

II.  EVALUATION  OF  METHOD  OF  MARAIS  AND  SHAW^®\ 

By  assuming  that:  (1)  stabilization  pond  contents  are 
completely  mixed,  (2)  decomposition  or  die-away  follow  first-order 
kinetics,  (3)  dally  average  values  of  flow  and  of  concentration  are 
sufficiently  accurate,  and  (4)  that  evaporation  can  be  Ignored, 
Marais  and  Shaw  derived  Equation  1.  (Also  derived  as  Equation  (A3f) 
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In  the  Appendix  of  thle  report) 

C  =  or  =  kT+1  (1) 

kT+1  C 

in  which 

=  concentration  or  number  in  Influent  flow 
C  =  concentration  or  number  in  pond  oont3nt8 
and  consequently  in  effluent  flow 
T  =  detention  time 

k  =  reaction  velocity  coefficient  for  a  first-order 
reaction 

This  equation  was  used  by  them  as  an  appropriate  mathematical  model 
to  describe  stabilization  pond  kinetics. 

They  collected  considerable  bacterial  count  data  over  periods 
of  months  from  several  ponds  in  Pretoria^  South  Africa  and  in 
Northern  Rhodesia.  The  bacteria  counts  Included  total  coliformsi 
E.  coll  and  F.  streptococci.  The  data  was  plotted  as  the  ratio  of 
the  initial  number  of  organisms  to  the  observed  number,  (.Ci/C) ,  at 
some  time  later  (detention  time)  as  the  ordinate  and  detention  time 
as  abscissa.  According  to  Equation  (1)  above,  the  data  should 
describe  a  straight  line  with  an  intercept  of  unity  and  a  slope 
equal  to  k.  There  appeared  to  be  quite  a  scatter  of  points;  however 
they  drew  the  straight  lines  and  values  for  k  were  taken  as  follows; 
Total  conform  k  =  2.13  day**^ 

E,  coll  2.14  >♦ 

F.  strep  2.82  ** 

For  purposes  of  design  k  was  taken  as  2.0.  The  design  equation  was 
then 
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No/N  =  2,0  T  +  1 
In  which  No  =  initial  count 

N  =  count  at  time  T 
T  -  time,  days 

Similarly  Marais  and  Shaw  collected  B,0.D.  (5  day-20°C}  data 
from  several  ponds  in  South  Africa,  N,  Rhodesia,  S,  Rhodesia  and 
from  several  pond  studies  in  the  United  States  Including  those  at 
MoJave^®\  Syracuae^^\  and  Fayette^^^,  The  foreign  data  was 
selected  from  regions  (or  seasonal  periods)  in  which  the  climatic 
conditions  were  similar  to  those  in  Southern  Africa.  It  was  in¬ 
tended  that  the  data  from  ponds  in  which  the  change  in  loading  rate 
and  sampling  schedules  were  Inadequate  would  not  be  included. 

However,  as  shown  below  some  of  the  data  from  the  United  States  was 
not  very  good.  The  data  was  taken  from  ponds  that  received  raw 
sewage,  settled  sewage  and  aqua  privy  effluents  (probably  septic 
tank  effluents).  The  depth  of  the  ponds  varied  from  two  to  ten  feet. 
Some  45  observations  were  plotted  according  to  Equation  (1)  above 
and  the  values  of  k  determined  to  be  0.23  day**^  for  Southern  and 
Central  Africa  and  0.17  day**^  for  the  U.  S,  data.  This  writer 
suggests  that  if  straight  lines  were  drawn  to  bracket  the  data,  the 
values  of  k  would  vary  from  about  0.06  to  0.4  day~^,  3y  arbitrarily 
excluding  a  couple  of  points  the  minimum  value  of  k  could  be  in¬ 
creased  to  0.1  day"^. 

One  of  the  major  objectives  of  this  present  paper  is  to 
evaluate  the  applicability  of  Equation  (1)  by  determining  the 
constancy  of  k.  Or,  it  is  suggested  that,  if  k  cannot  be  shown 
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to  be  reasonably  eonstantf  tben  the  validity  of  Equation  (1)  as  a 
mathematical  model  is  to  be  reconsidered,  Aoeordingly,  considerable 
data  were  collected  (it  is  believed  almost  all  available  and  applic¬ 
able)  from  various  pond  studies  in  the  United  States,  The  data  is 
tabulated  In  Tables  1  through  9  and  plotted  according  to  Equation 
(1)  on  Figures  1  through  9,  The  discussion  follows. 

Data  collected  at  Fayette,  Missourl^^^  are  presented  in  Table  1 
together  with  calculated  values  of  k  based  on  Equation  (1),  This 
data  Is  also  plotted  on  Figure  1.  Inspection  of  these  results  re¬ 
veal  that  overall  rate  constant  kt 

1,  varied  from  0,044  to  0,759, 

2,  was  relatively  insensitive  to  temperature  variations, 

3,  varied  approximately  inversely  with  detention  time, 

(  S  ) 

Data  collected  at  Farmvllle,  Virginia  are  presented  In 
Table  II  together  with  calculated  values  of  k  based  on  Equation  (1), 
This  data  is  plotted  on  Figure  2,  Inspection  of  these  results  re¬ 
veal  that  the  overall  rate  constant  k  for  any  single  ponds 

1,  varied  from  0,00  to  0,135, 

2,  increased  from  Itiarch  to  June  and  decreased  during  July 
and  August,  During  September  and  October  the  average  k 
In  ponds  A  and  C  compared  favorably  with  the  March  to 
August  average. 

Average  k  during  Average  k  during 
March  to  Aug, _  gtPJfes.  lft-22Li _ 

POND  A  0,033  0.039 

POND  C  0,078  0,092 
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TABLE  I  BOD  REMOVAL  AT  FAYETTE,  MIS: 


PERIOD 

INFLUENT 

BOD 

r 

EFFLUENT 

■  BOD  C  -  ^ 

HG/L 

Ci 

MG/L 

CELL  1  ^ 

CELL  2  ® 

I 

CELL  3  ^ 

CELL  4  ^ 

CELL  5^ 

1 _ 

CEL 

1957 

May 

254 

35 

31 

38 

34 

59 

0.0' 

Jun 

514 

32 

36 

50 

54 

39 

1  o,at 

Jul 

279 

44 

46 

37 

45 

4l 

o.oi 

Aug 

266 

43 

37 

53 

57 

55 

o.oi 

Sept 

508 

45 

56 

57 

52 

51 

0.01 

Oct 

280 

3c. 

47 

55 

51 

0.0' 

Hov 

266 

33 

53 

44 

29 

38 

0.0' 

Dec 

310 

34 

57  1 

43 

37 

1  4o 

0.0« 

1958 

Jan 

270 

27 

32 

75 

71 

42 

0.1: 

Feb 

223 

28 

40 

50 

63 

63 

o.o< 

Mar 

252 

51 

34 

48 

57 

66 

o.o< 

Apr 

195 

31 

32 

42 

42 

4o 

0.0' 

Hay 

252 

28 

33 

37 

40 

51 

o.r 

a  CjII  1  -  depth  2.5  ft,  area  0.7^  acres,  leading  L  -  20.5  Ib/acre/day,  average  detenl 

b  Cell  2  -  depth  2.3  ft,  area  0.75  acres,  loading  L  ■  40.5  Ib/acre/day,  average  detonl 

I  c  Cell  5  -  depth  2.5  ft,  area  0.75  acres,  loading  L  -  CO.F  Ib/acrc/day,  average  detenl 

I  d  Cell  -  depth  2.3  ft,  area  0.75  acres,  loading  L  »  8l.l  Ib/acro/day,  average  detenl 

e  Cell  5  -  depth  2.3  ft,  area  0.75  acres,  loading  L  -101.5  Ib/acre/day,  average  detenl 


CELL  1^ 

CELL  2® 

CELL  3*^ 

CELL  4  ' 

CELL  S'J 

1958 

Jul 

253 

25 

27 

42 

22 

35 

1  0.7' 

Aug 

255 

39 

21 

38 

26 

37 

o.3i 

1  Sept 

285 

51 

18 

39 

38 

44 

0.3; 

Oct 

285 

61 

21 

43 

53 

36 

0.2; 

Hov 

289 

51 

22  1 

46 

37  ; 

55  I 

0.3! 

Dec 

286 

58 

25  ! 

1  51 

43  1 

38 

0.2‘ 

:  1959 

Jan 

267 

69 

25 

55 

48 

51 

0.1' 

1  Fob 

272 

77 

4o 

•  67 

52 

54 

o.r 

Mar 

245 

55 

36 

66 

51 

!  56 

0.2; 

Apr 

247 

47 

'  1 

30 

49 

52 

52 

0.? 

f  Cell  1  -  depth  2.5  ft,  areaa.75  acres,  loading  L  ai20  Ib/acre/day,  average  detention 
g  Cell  2  -  depth  2.5  ft,  area  0.75  acres,  loading  L  a  receives  effluent  from  Cell  5,  a 
b  Cell  5  -  depth  2.5  ft,  area  0.75  acres,  loading  L  ^00  Ib/acre/day,  average  detoutloi 
i  Cell  4  -  depth  3.0  ft,  area  1.00  acres,  loading  I  ■  60  Ib/aere/day,  average  detent  la 
J  Cell  5  -  depth  2.5  ft,  area  0.75  acres,  loading  I  ■  60  Ih/acre/day,  average  detentim 
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TABLE  I  BOD  REMOVAL  AT  FAYETTE ,  MISSOURI 


EFFLUENT  BOD  C  -  MG/L 

K  -  days"' 

1 

CELL  2  ® 

— 

CELL  3 

CELL  4  ^ 

CELL  5  ^ 

CELL  1 

CELL  2 

— 

CELL  3 

CELL  4 

CELL  5 

31 

38 

34 

39 

0.077 

0.177 

0.210 

0.319 

0.338 

36 

50 

54 

39 

0.087 

0.153 

0.157 

0.191 

0.349 

46 

37 

45 

41 

0.048 

0.090 

0.177 

0.187 

0.261 

37 

53 

57 

53 

0.044 

0.105 

0.102 

0.124 

0.170 

56 

57 

52 

51 

0.055 

0.085 

0.125 

0.186 

0.239 

47 

55 

51 

50 

0.075 

0.109 

0.135 

0.198 

0.253 

33 

44 

29 

38 

0.079 

0.159 

0.170 

0.368 

0.339 

37 

43 

37 

4o 

0.093 

0.170 

0.214 

0.339 

0.388 

32 

75 

71 

42 

0.111 

0.182 

0.096 

0.138 

0.333 

40 

50 

63 

63 

0.090 

O.U8 

0,134 

0.130 

0.163 

34 

48 

57 

66 

0.093 

0.167 

0.166 

0.178 

0.183 

32 

42 

42 

40 

0.076 

0.147 

0.157 

0.209 

0.279 

33 

37 

40 

51 

0.117 

0.195 

0.256 

0.310 

0.289 

.7?  acres,  loading  L  -  20.5  Ib/acre/day, 
.75  acres,  loading  L  ■  40.  b  Ib/acre/day, 
.75  acjres,  loading  L  ■  CO.F  Ib/acro/day, 
.75  acres,  loading  L  »  fil.l  Ib/acro/day, 
.75  acres,  loading  L  -101.5  Ib/acre/day, 


average  detent! cn  time  T 
average  detcntioo  tine  T 
average  detention  time  T 
average  detention  tine  T 
averSLge  detention  tine  T 


>  67  days 
-  44  days 

■  29  days 
«  22  days 

■  17  days 


CELL  2^ 


CELL  3  CELL  4 


CELL  5'^ 


27 

21 

18 

21 

22 

25 

25 

40 

% 

50 


42 

38 

39 

46 

51 

55 

67 

66 

49 


22 

35 

0.759 

1  0.296 

0.135 

0.215 

26 

37 

0.382 

0.336 

0.113 

0.203 

38 

44 

0.516 

0.371 

0.083 

0.189 

33 

36 

0.255 

0.331 

0.098 

0.239 

37 

35 

0.522 

— 

0.511 

0.087 

0.250 

43 

58 

0.271 

— 

0.271 

0.072 

0.225 

48 

51 

0.198 

w 

0.227 

0.058 

0.146 

52 

54 

0.175 

•  • 

0.173 

0.054 

0.139 

51 

56 

0.236 

0.158 

0.048 

0.115 

52 

52 

0.2^ 

1  0.238 

0.048 

0.129 

'5  acres,  loading  L  ^20  Ib/acre/day,  average  detention  T  -  14.5  days 

1.75  acres,  loading  L  .  receives  efl’luent  from  Cell  5,  average  detention  T  >  54  days 

1.75  acres,  loading  l  ^00  Ib/acre/day,  average  dotoution  T  -  17  days 

.00  acres,  loading  I  -  60  Ib/acre/day,  average  detention  T  -  78  days 

1.75  acres,  loading  I  -  60  Ib/acre/day,  average  detention  T  -  29  da^ 
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FIGURE  I  BOD  REMOVAL  AT  FAYETTE,  MISSOURI 
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TABLE  m  BOD  REMOVAL  IN  WISCONSIN  OXIDATION  PONDS 
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FIGURE  5  BOO  REMOVM.  m  PILOT  PLANT  OXIDATION  PONDS  AT  SYRACUSE,  NEW  YORK 
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FIGURE  e  BOO  REMOVAL  AT  MOJAVE,  CALIFORNIA 
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TABLE  m  BOD  REDUCTION  AT  RICHMOND,  CALIFORNIA 


PERIOD 

INFLUENT 

BOD 

Ci 

MG/L 

EFFLUENT 

BOO 

C 

M6/L 

DETENTION 

TIME 

T 

DAYS 

DEPTH 

0 

INCHES 

LOADING 

L 

LB/ACRE/OAY 

K-  DAYS"' 

PGR)  I 

Ihy  1954 

226 

43 

3 

8 

136 

1.42 

Jul  1954 

150 

12 

3 

8 

93 

3*63 

Au«  1954 

115 

25 

3 

8 

71 

1*35 

npy  1956 

116 

29 

30 

36 

32 

o.ib 

Jul  1956 

141 

32 

20 

36 

52 

0.17 

Aug  1956 

117 

35 

10 

36 

106 

0*23 

Sapt  1956 

171 

58 

5 

36 

232 

0.39 

PORD  II 

JUl  1954 

200 

A6 

4 

12 

115 

0.64 

Aug  1954 

173 

19 
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12 

153 

2.70 

Aug  1954 

175 
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2 

12 

231 

2*63 

Aug  1954 

171 
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12 

440 

1.04 

PdD  111 
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36 
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14 

114 

1.14 

Aug  1955 

153 

46 
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0.58 
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25 
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0.61 
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22 
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1.42 
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FIGURE  7  BOO  REMCWL  AT  RICHMOND,  CALIFORNIA 
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TABLE  H  COLIFORM  REMOVALS  IN  AUSTRALIA  AND  CALIFORNIA 


PERIOD 

INFLUENT 

E.  COLI./  ML 

c: 

EFFLUENT 

E.  COLI./ML 

C 

DETENTION 

TIME 

T 

DAYS 

K'DAYS"' 

Auatralla  *'  Murtcaioi  -  Aerobic  Pond 

Summer 

35 » 000 

5,874 

10.5 

0.8 

Vinter 

98,750 

860 

20.0 

5.7 

Australia  •  115B  *■  Aerobic  Pond 

Sumner 

51,000 

5.100 

17.5 

0.9 

Vinter 

89,100 

5,600 

37.0 

0.6 

Auatralla  *■  Murtcala  -  Anaerob 

Le  Pond 

Sumner 

125,000 

55,000 

3.5 

0.7 

Vinter 

162,000 

55,000 

3.5 

0«6 

Auatralla  *•  115B  "  Anaerobic  Pond 

Suonar 

1.05x10^ 

143,000 

3.5 

1.8 

Vinter 

162,000 

89,000 

7.0 

0.1 

California  •  Vew  Pond 

1955-56 

Bov  26-Jan  16 

61,000 

30,000 

7 

0.15 

Jan  16-Fab  20 

72,000 

15,000 

30 

0.13 

Peb  pOHfer  Vf 

55,000 

8,300 

43 

o.oe 

far  13-fcr  27 

67,000 

13,000 

30 

0.14 

■ar  ih’-kpr  2* 

45,000 

32,000 

9 

0.05 

Apr  24-lley  15 

150,000 

13,000 

17 

0.62 

California  -  Vest  Lagoon 

Sept  27-Oct  26 

80,000 

61,000 

1.2 

0.26 

Oct  2e<an  16 

40,000 

24,000 

6.9 

0.10 

Jan  16Hiar  27 

60,000 

13,000 

1.3 

2.79 

■ar  27-Apr  2* 

45,000 

55,000 

2.4 

Apr  2dHlay  15 

150,000 

26,000 

1.8 

2.65 

24 


The  k  dropped  by  60  to  64  per  cent  from  September 
to  October  In  ponds  A  and  C, 

3.  varied  approximately  Inversely  with  detention  time. 

Data  collected  during  the  studies  In  Wisconsin  oxidation 

(  O  \ 

ponds  are  presented  In  Table  III.  Figure  3  Is  a  plot  of  this 
data « 

The  overall  rate  constant  ki 

1.  varied  from  0.011  to  1.359. 

2.  was  higher  during  the  summer. 

3.  did  not  vary  Inversely  with  detention  time  although  a 
comparison  between  the  New  Auburn  and  Spooner  ponds 
during  August  revealed  that  k  decreased  with  Increasing 
detention  time. 

_  (7) 

Data  collected  In  pilot  plant  oxidation  ponds  at  Syracuse 
are  presented  in  Table  IV  and  Table  V,  Figures  4  and  9  are  plots 
of  the  data.  The  data  shows  considerable  variation  and  conclusions 
are  difficult  to  draw.  Part  of  the  variation  can  be  attributed  to 
the  small  number  of  samples  which  were  analyzed.  Only  one  or  fewer 
samples  were  collected  each  week.  Table  V  describes  conditions  at 
Syracuse  during  the  last  five  weeks  of  the  summer  when  the  loadings 
were  increased  weekly  from  104  Ib.BOD/acre/day  to  1246  Ib.BOD/aore/ 
day.  The  variation  in  k  was  extreme,  values  varying  from  0.026  to 
6.26.  It  is  evident  that  loadings  in  some  basins  were  changed 
before  equilibrium  conditions  had  been  attained.  Failure  to  con¬ 
tinue  a  constant  loading  for  a  sufficient  period  resulted  in  samples 
being  collected  which  represented  the  tail-end  of  the  previous  weeks 
loading. 
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Data  collected  at  Mojave,  Callfornla^^^  are  found  in  Table  VI 

and  plotted  on  Figure  6|  and  data  collected  at  Richmond,  Cali* 

(3) 

fornia  in  Table  VII  and  Figure  7.  At  Mojave,  the  overall  rate 
constant  k  appears  to  be  Inversely  proportional  to  detention  time 
with  values  of  k  ranging  from  0.03  to  0.39,  At  Richmond  k  varied 
from  0.1  to  3.83  and  showed  a  time  relationship  only  for  detention 
times  greater  than  ten  days. 

The  available  data  on  the  removal  of  coliform  organisms  in 

oxidation  ponds  is  more  limited  than  BOD  removal  data.  Data  from 

(4) 

Fayette,  Missouri  is  found  in  Table  VIII  and  plotted  on  Figure  8. 
Data  from  Australia^ and  California^ are  found  in  Table  IX 
and  plotted  on  Figure  9.  These  data  showed  extreme  variation  with 
rate  constants  varying  from  0.05  to  1,449,  There  was  no  correlation 
of  coliform  reduction  and  detention  time. 

Thus  it  appears  upon  more  extensive  evaluation  that  Equation 
(1)  does  not  serve  very  well  as  a  mathematical  model  to  describe 
the  decrease  of  a  waste  constituent  with  time.  The  data  from  the 
several  sources  reveal  generally  that  the  value  of  k  varies  widely 
for  no  readily  apparent  reason,  is  not  dependent  upon  temperatures 
above  freezing,  and  is  approximately  Inversely  proportional  to 
detention  time.  An  attempt  is  made  in  the  succeeding  Section  to 
develop  another  relationship, 

III.  DISCUSSION  OF  DISCREPANCIES  BETWEEN  PROPOSED  THEORY 
AND  EXPERIENCE. 

It  is  apparent  from  the  data  analysis  in  section  II  of  this 
report  that  the  proposed  means  of  evaluating  pond  performance  has 
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certain  deficiencies.  It  is  the  purpose  of  this  section  of  the 
report  to  suggest  reasons  for  the  apparent  deficiencies  and  to 
propose  another  means  of  evaluation  of  pond  performance. 

The  commonly  used  detention  periods  in  stabilization  ponds  is 
several  times  the  period  normally  required  for  essentially  complete 
biochemical  oxidation  of  the  putresclble  organic  matter  in  sewage. 
Because  of  the  biological  system  in  a  stabilization  pond,  the 
organic  matter  in  the  sewage  is  decomposed  by  bacterial  action  and 
the  nutrients  released  are  converted  in  part  to  the  plankton, 
primarily  algae,  and  a  biological  cycle  is  established  —  as  the 
algae  die,  the  bacteria  decompose  the  dead  cells  releasing  nutrients 
for  more  algae,  etc.  The  continuous  source  of  energy  is  the  sun, 
and  the  continuous  sewage  additions  provides  the  replacement  of 
dissipated  or  discharged  nutrients.  The  oxygen  released  by  the 
photosynthetlo  processes  of  the  algae  are  utilized  by  the  bacteria 
to  decompose  aerobically  the  organic  matter  present,  A  state  of 
equilibrium  may  be  reached  if  constant  conditions  of  sewage  feed, 
sunlight,  temperature,  etc.  persist  for  a  period  of  time.  The 
effluent  from  such  a  pond  would  contain  organic  matter,  to  some 
measure  living  on  dead  plankton  cells,  and  if  the  ordinary  means 
of  BOD  measurement  is  used,  the  decay  of  these  cells  in  a  bottle 
in  the  dark  would  yield  an  oxygen  demand.  Depending  upon  the  form 
or  state  of  the  receiving  water  course  or  land  surface  the  effect 
of  the  living  algae  cells,  which  are  measured  as  BOD,  may  be  sig¬ 
nificantly  different  than  an  equivalent  amount  of  BOD  in  the  form 
of  say  sewage.  However,  the  BOD  is  then  calculated.  (It  is  not 
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tv,-,  tM<8  report  to  do  so,  but  it  might  be  Interesting 

to  try  to  determine  what  portion  of  the  BOD  of  a  pond  effluent  is 
in  the  form  of  plankton  cells  and  what  portion  is  undecomposed 
sewage,) 

It  is  suggested  that  changes  in  the  theoretical  detention  times 
of  a  stabilization  pond  will  not  result  in  a  corresponding  change 
in  the  BCX)  of  the  effluent.  Thus,  evaluation  of  the  rate  of  sewage 
decomposition  by  measurement  of  the  effluent  BCD  and  relating  those 
values  to  the  Influent  BOD  may  not  be  as  straightforward  a  procedure 
as  suggested  by  Marais  and  Shaw^^\  For  example,  if  say  20  days  is 
sufficient  time  for  a  pond  to  reach  a  certain  value  of  effluent  BCD, 
Increase  of  the  detention  time  to  say  40  days  might  not, as  suggested 
above,  result  in  a  change  in  effluent  BQD»  However,  the  value  of  k 
will  be  changed  by  a  factor  of  approximately  two.  Or,  after  a 
certain  minimum  detention  time  in  a  pond  (value  not  yet  determined) 
increase  in  detention  times  would  result  in  a  decrease  in  k.  This 
is  Illustrated  in  Figures  1-9;  if  C  reaches  an  equilibrium  values, 
Ci/C  is  constant  and  the  slope  k  changes  with  detention  time  T. 

With  this  concept  in  mind,  the  values  of  k  determined  as  the 
slopes,  are  plotted  vs,  detention  time  T  in  Figures  10  and  11,  It 
is  seen,  most  clearly  on  Figure  10,  that  the  data  seem  to  describe 
a  hyperbola  -  or  k  T  might  be  considered  a  constant.  From  Equation 

(  1 )  ,  P 

^1 

C  =  Z -  ,  (1) 

kT+1 

it  is  seen  that  if  kT  is  a  constant,  C  is  a  constant  fraction  of  . 
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OeTENTON  TIME  -  OWS 

FIGURE  10  BOO  REMOVHl.  RATE  CONSTANT  VERSUS  DETENTION  TIME  AT  FAYETTE,  MISSOURI  AND  FARMVILLE,  VIRGINIA 


31 


FIQURC  It  800  REMOVAL  RATE  CONSTANT  VERSUS  DETENTION  TIME 


Or,  Irrespective  of  detention  tlae  the  sane  fraction  removal  is 
effected , 

To  check  this  assumption,  values  of  BOD  removal  (R)  in  terms 

of  mg/l/day  as  shown  on  Figures  1^  -  1^.  On  Figure  12,  the  plot 

(4)  - 

of  the  Fayettej  Missouri  data  ,  the  average  ratio  of  R/L  equals 

approximately  0.85  and  the  range  of  detention  time,  from  Figure  1 

is  10  to  160  days.  This  removal  efficiency  was  as  effective  at  the 

maximum  loading  rate  of  about  20  mg/l/day  as  it  was  at  the  lower 

loading  rates.  Figure  13  is  a  plot  of  the  data  from  the  Texas 

studies^  the  ratio  of  R/L  has  an  average  value  of  about  0.9  and 

the  maximum  loading  is  approximately  90  mg/l/day.  On  Figure  14, 

the  data  from  the  Farmville,  Virginia  studies  ,  the  average  value 

of  R/L  is  about  0.85.  On  Figure  15,  the  data  is  from  the  Wisconsin 
/  q  \  —  — 

studies  ,  the  ratio  of  R/L  is  about  0,80  except  at  the  highest 
loading  when  the  ratio  is  1.0.  This  last  point  doesn't  seem  to 
have  any  significance.  Thus  we  see  that  ratio  of  R/L  lies  in  a 
rather  narrow  range  of  0.80  to  0.90,  More  significantly,  perhaps, 
is  that  for  a  given  pond  the  ratio  is  very  nearly  constant,  l.e, 
a  straight  line  seems  to  represent  the  data  very  well.  A  signifi¬ 
cant  increase  in  loading  rate,  over  those  observed  from  these  data, 
might  result  in  a  break  in  the  line  relating  L  to  R  on  Figures 
12-15.  That  is,  if  the  efficiency  of  removal  drops  the  slope  of 
the  line  would  decrease. 
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FIGURE  12  BOO  REMOVAL  VERSUS  BOO  LQMMNG  AT  FATETTE,  MISSOURI 
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FIQURE  »  BOD  REMOVAL  VERSUS  BOD  LOAOINe  AT  WISCONSIN 


IV.  POSSIBLE  AREAS  FOR  FURTHER  STUDY 


In  section  II  of  this  report  It  was  shown  that  the  proposed 
relationships between  the  ratio  of  incoming  to  outgoing  BOD  con¬ 
centrations  or  bacteria  numbers  as  a  linear  function  of  detention 
time  was  not  a  good  representation  of  the  accumulated  data.  In 
section  III  the  primary  reason  for  such  deviations  as  noted  In 
section  II  was  proposed;  namely  that  BCX)  reduction,  and  possibly 
reduction  In  bacteria  numbers,  were  essentially  completed  In  a 
period  of  time  much  less  than  the  detention  periods  commonly  allowed 
In  ckldation  ponds.  However,  there  must  be  an  upper  limit  of  load¬ 
ing  rate,  or  reduction  In  detention  time,  at  which  the  efficiency 
of  removal  Is  reduced.  This  Is  Illustrated  by  Figure  16. 

Removal 

mass/unit  volume/unit  time 


Loading,  mass/unlt  volume/unit  time 
FIGURE  16 

Relationship  between  Loading  and  Removals  In  an  Oxidation  Pond 
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The  rising  part  of  the  curve  in  Figure  16  is  illustrated  by  real 
data  on  Figures  12-19.  The  data  'available  for  this  study  did  not 
reach  to  the  region  where  the  curve  flattens  out.  That  is,  the 
critical  loading  value  had  not  been  reached. 

If  one  is  primarily  interested  in  removal  efficiencies,  opera¬ 
tion  of  the  pond  in  any  region  up  to  the  critical  loading  value 
would  yield  an  effluent  with  a  concentration  of  say  0.1  to  o«2 
the  concentration  of  the  influent,  l.e.  C^/C  =  5  to  10,  or  per¬ 
centage  removals  are  80  to  90%.  In  this  loading  region,  l.e.  the 
rising  part  of  the  curve,  the  range  of  removal  efficiencies  is  due, 
not  to  differences  in  detention  time,  but  rather  to  differences  in 
conditions  such  as  temperature,  pH,  presence  or  absence  of  essen¬ 
tial  nutrients,  sunlight  intensity  and  duration,  degree  of  cloudi¬ 
ness  etc.  Some  of  these  above  listed  conditions  vary  with  the 
season.  Even  so,  variation  in  such  conditions  may  have  a  reasonably 
small  effect  cn  the  slope  of  the  rising  part  of  the  curve  on  Figure 
16,  but  might  shift  the  break-off  point  or  critical  loading  value 
rather  substantially.  Determination  of  which  effect  would  be  most 
likely,  i.e.  change  in  slope  or  location  of  critical  loading  value, 
would  be  an  interesting  topic  for  further  study. 

The  suggestion  of  single  pond  removals  Illustrated  by  Figure  16 
could  well  apply  to  ponds  in  a  series.  That  is,  if  n  second  pond  could 
remove  say  80-90%  of  the  influent  BOD  at  any  loading  leas  than  the 
critical  loading  value,  dividing  it  up  into  several  ponds  in  series 
would  lead  to  significant  overall  increase  in  efficiency,  Marais 

_  (5) 

and  Shaw  have  illustrated  this  concept  and  many  pond  designers 
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have  recommended  it. 

Of  course  the  discussion  to  this  point  has  not  included  con¬ 
sideration  of  the  aerobic  or  anaerobic  conditions  of  the  pond, 
Oswald^ suggests  that  most  operating  oxidation  ponds  are  an¬ 
aerobic  or  facultative.  The  latter  is  anaerobic  near  the  bottom 
and  aerobic  nehr  the  top.  In  the  anaerobic  state  'good  methane** 
fermentation  may  block  reduction  of  sulfur  compounds;  thus  a  state 
of  anaeroblases  doe'  not  necessarily  result  in  an  odor  nuisancSf 
particularly  if  the  top  layers  of  the  pond  are  aerobic.  Thus, 
one  cannot  dismiss  the  concept  of  the  critical  loading  value  as  a 
design  criteria  on  the  basis  that  anaerobic  conditions  might  exist. 
In  fact  during  several  hours  per  year,  if  not  several  hours  per  day, 
all  oxidation  ponds  might  be  anaerobic  through  the  entire  depth, 
without  creating  odor  problems. 

One  additional  concern  in  loading  must  be  discussed,  that  is 
the  condition  of  icing.  In  the  northerly  regions  of  the  United 
States  and  of  course  other  similar  climatic  regions,  an  oxidation 
pond  will  be  frozen  over  for  several  days  to  several  months  each 
year.  During  this  time  biological  decomposition  has  been  slow 
and  organic  matter  has  accumulated.  However,  due  to  sedimentation 
overall  removals  may  be  good.  Upon  warming  biolcpioal  decomposition 
begins  and  is  probably  anaerobic.  During  this  spring  breakup  period 
free  oxygen  made  available  by  photosynthesis  or  aeration  is  in¬ 
adequate  to  maintain  aerobic  conditions.  The  pond  or  aeration  may 
emanate  odors.  The  greater  the  accumulation  of  organic  matter,  and 
this  is  of  course  related  to  loading  rates,  the  longer  the  odor  con- 
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dltlon  will  persist.  This  too  might  be  an  interesting  area  for 
study . 

Thus  I  it  appears  that  as  operating  and  experimental  evidence 
is  accumulated,  bases  for  design  may  be  reconsidered,  and  this  writer 
suggests  that  the  method  of  expressing  results  suggested  by  Figure 
16,  that  is  removals  vs.  loadings  both  in  terms  of  naas/unlt  volume/ 
unit  time,  may  prove  to  be  of  value. 
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APPENDIX 


If  a  naterlal  balance  is  aade  over  the  oxidation  pond  the  word 
equation  yields  Equation  (A- 1), 

Inflow  -  Outflow  +  Production  =  Aeeuaulatlon  (A~l> 

If 

Qj  =  Inflow  voluaetric  rate,  (L^T**^) 

Qo  =  effluent  voluaetric  rate,  (L3T"1) 

V  -  pond  voluae  L 

Cji^  =  Inflow  concentration  (ML*’3} 

Cj,  =  concentration  at  any  point  in  the  pond  (ML'*3) 

C  =  effluent  concentration  (ML'S), 

Equation  (1)  can  be  written 

7  ^7 

QiCi  -  QoC  -  /^RdV  =  -57  ^  Cl  dv  (a-2) 

Note  that  the  production  tera  has  a  negative  sign  to  indicate  de- 
orease  in  waste  constituent  concentration*  "K"  represents  the  rate 
at  which  waste  constituent  oonoentratlon  is  being  reaoved  by  physi¬ 
cal,  cheaioal  or  biological  neans  and  has  the  units  of  * 

Since  the  reaction  rate  nay  depend  upon  the  concentration  Cl  and 
since  Cl  nay  vary  fron  point  to  point  in  the  pond,  K  aust  be  written 
inside  the  integral*  For  tho  sane  reason  the  aeeuaulatlon  tera 
of  Equation  (A>3),  which  represents  the  total  change  in  mass  with  tine, 
aust  be  written  as  the  derivative  of  an  Integral*  In  fact,  all  the 
quantities  in  Equation <A-2)  any  vary  with  time. 

Equation  (A-2)  nay  be  modified  for  various  conditions*  Several 
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different  cases  are  discussed  below 


Case  1  -  Lagoon  Contents  Conpletely  Mixed, 

The  concentration  at  any  point  in  the  pond  is  equal  to 
the  effluent  concentration,  l,e,  =  C.  Equation  (A-2) 

can  then  be  written  as 

QjC  -  Q^C  -  KV  =  (A-3) 

Note  that  the  integrals  do  not  appear  in  Equation  (A-3) 
because  the  effluent  concentration  C  is  a  function  only 
of  the  total  volume  V  and  not  a  function  of  position  with 
pond , 

Case  1-a  -  Complete  Mixing  and  Constant  Voluiaetrlo  Flow  Rate. 

If  =  Qq  Equation  (A-3)  simplifies  to  Equation  (A-3a) 

(i(C,-C)  -  KV  =  (A-3a) 

dt 

Case  1-b  -  Complete  Mixing,  Constant  Mass  in  Flow  Rate,  Constant 
Pond  Volume  and  Steady  State. 

In  this  case  dC/dt  =  0  and  Equation  (A-3)  can  be  written 
Q(Ci-C)  -  KV  =  0,  (A-3b) 

If  the  order  of  the  reaction  is  known  the  relationship  between 
time  and  effluent  concentration  can  be  derived.  For  example, 
in  the  case  of  a  concentration  reduction  in  accordance  with 
a  first-order  reaction, 

K  =  kCi 

so.  Equation  (A-3b)  can  be  written 

Q(Cj-C)  -  kCn^  =  0  (A-3c) 
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to  yield  Equation  (A-3e). 

By  dividing  though  by  V,  one  obtains 

S  (C^  -  C)  -  kC  =  0.  (A-3d) 

Since  V/Q  -  detention  time  T,  Equation  (A-3d)  may  be 
written 


(Cj  -  C)  =  kCT  (A-3e) 

or 

^1 

C  =  +  1  (A3£) 

which  Is  also  Equation  (1)  In  the  body  of  this  report 
and  the  working  equation  derived  similarly  by  Marias 
and  Shaw^  , 

Case  2  -  Steady  State  Conditions,  Imperfect  Mixing, 

If  Qj  =  Qq,  Equation  (A-2)  can  be  wrv'itten 

V 

Q<Ci  -  Cq)  -  /  kdV  =  0,  (A-4) 

0 

If  only  overall  conditions  are  of  Interest  the  Integral 
can  be  replaced  by  RV  to  yield, 

Q<Cj^  -  Cq)  -  RV  =  0  (A-4a) 


or 


(A-4b) 


R  is  in  the  removal  effected  by  the  pond  in  terms  of 
mass/unit  time/unit  volume.  Relating  R  to  L,  the  loading 
on  the  pond,  (L  =  QC^/V)  may  allow  evaluation  of  the  pond 
without  knowledge  of  the  degree  of  mixing,  extent  of 
sedimentation,  or  reaction  order.  The  results  of  such 
evaluation  of  R  vs.  L  are  shown  on  Figs,  12-19  in  the  body 
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of  this  report 


Incidentally  Equation  (A-3b)  could  be  obtained  from 
Equation  (A-da)  by  assuming  a  zero-order  reaction  and  simply 
designating  R  as  X. 
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